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The effect of antiferromagnetic interchain coupling in alternating spin (1,1/2) chains is studied by 
mean of a spin wave theory and density matrix renormalization group (DMRG). In particular, two 
limiting cases are investigated, the two-leg ladder and its two dimensional (2D) generalization. Re- 
sults of the ground state properties like energy, spin gap, magnetizations, and correlation functions 
are reported for the whole range of the interchain coupling J± . For the 2D case the spin wave results 
predict a smooth dimensional crossover from ID to 2D keeping the ground state always ordered. For 
the ladder system, the DMRG results show that any J± > drives the system to a gapped ground 
state. Furthermore the behaviour of the correlation functions closely resemble the uniform spin-1/2 
ladder. For J± lower than 0.3, however, the gap behaves quadratically as A ~ Finally, it is 

argued that the behaviour of the spin gap for an arbitrary number of mixed coupled spin chains is 
analogous to that of the uniform spin-1/2 chains. 

PACS numbers: 75.10-b,75.10Jm,75.50Ee,75.50.Gg,76.50.+g 



I. INTRODUCTION 

The recent synthesization of quasi-onedimensional 
bimetallic magnets ]l],§, with each unit cell containing 
two kind of different spins, has spurred a growing interest 
in the low-temperature properties of quantum ferrimag- 
netic chains. After intensive studies it was possi- 

ble to describe properly the thermodynamic properties 
observed in experiments by including the ferro and anti- 
ferromagnetic (AF) features that make this systems spe- 
cially attractive. These results together with the possible 
experimental realization of new mixed spin compounds 
motivated the study of a variety of AF systems with 
different spin composition, all kind of interactions, and 
topologies (h]-[k|. From the experimental point of view, 
the ferrimagnetic compounds found are all composed of 
weakly coupled alternating spin chains Q| (Figj^(b) rep- 
resents the ladder case of such a topology). In other fam- 
ily of compounds like M nC 'u{pba) (HzO)^. .2 H2O, with 
pba=l,3-propylenebis(oxamato), the ground state is a 
non-magnetic one, although it is composed of weakly 
coupled alternating spin (Mn II Cu 11 ) chains |l| (Figj^(a) 
and (c)). In fact, using the Heisenberg Hamiltonian for 
these interacting spin systems and applying Lieb-Mattis 
theorem [[id} , it is easy to prove that in case (a) and (c) 
the ground state is AF with total spin Stot — whereas 
in case (b) it is ferrimagnetic with Stot 7^ 0. Nonetheless, 
in the 2D case we will see that the rotational spin sym- 
metry can be broken in the thermodynamic limit. The 
aforementioned considerations motivated us to focus on 
the ground state properties of an alternating spin ladder 
of type (a) and a 2D array of type (c) with the following 
Hamiltonian, 

H = J S i- s J + ^ E (Si-Sj+Si.Sj), (1) 

<i,J>|| <i,J>x 



where Si(si) represents a spin-1 (spin-l/2)and < i,j >|| 
(< i,j >±) denotes nearest neighbors along horizontal 
(vertical) direction. To our knowledge, this is the first 
systematic study of antiferromagnetically coupled alter- 
nating spin chains. Only recently, this issue has been ad- 




FIG. 1. Three kind of alternating spin systems: (a) and (c) 
have a singlet ground state, (b) has a ferrimagnetic one. 



dressed but for ferrimagnetic ladders (Fig.|T](b) ) M,ff6[. 
Instead, in the alternating spin ladder we are interested 
in, (in what follows ASL), the global properties of the 
system change completely once the chains are coupled. 
Unlike uniform spin ladders, little is known about the ef- 
fect of AF interchain coupling in ASL. For instance, in 
a uniform spin-1/2 ladder (USL), it is well known that 
an infinitesimal J± is required to send the chains off- 
criticality -with an excitation gap growing linearly with 
J_i_- independently of its sign [pTyiq] , Besides, it has 
been pointed out that the absence of broken symmetry 
in the chains is indeed the crucial feature for such a be- 
haviour [jl7|]. On the other hand, in the spin-1 case, it 
has been argued that a interchain coupling J± makes the 
gap (already existing in the chains) to decrease with Jj_, 
and again no matter the sign of J± For the mixed 
case, we are concerned, the situation is rather different as 
the ferrimagnetic chains are neither gapped nor critical, 
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but they are gapless as a result of a broken symmetry 
ground state. Hence, one could intuitively say that the 
effect of interchain coupling should be different to the 
both uniform cases mentioned above and a finite amount 
of J± could be needed to melt order. We shall sec, how- 
ever, that the low dimensionality of the problem (ID) 
is crucial and that any Jj_ > drives the system to a 
gapped ground state. The study of the spin correlation 
functions confirms this picture, showing a strong resem- 
blance to USL. Contrary to the ASL, in the 2D case we 
shall see that the dimensional crossover from ID is com- 
pletely smooth and keeps the ground state ordered for 
the whole range of J±. 

The work is organized as follows: in section II we study 
the dimensional crossover from ID to 2D within spin 
wave theory. In section III we give some details of the 
DMRG method implemented and discuss the results for 
the ladder case. In section IV we summarize the results 
and give the concluding remarks. 



II. SPIN WAVE THEORY 

Though approximate, it has been shown that the spin 
wave (SW) series converges correctly to the DMRG re- 
sults in ferrimagnetic chains ||. Hence SW appears to 
be, in principle, a good starting point technique to study 
the crossover from ID to 2D where SW is even more reli- 
able. It is also known, however, that for a ladder system 
SW breaks down |L2j (see below). The problem requires 
a unit cell composed of two spins 5 = 1 and two spins 
5=1/2 (see Figj^(c)). At the classical level, it is as- 
sumed a Neel order and quantum fluctuations are incor- 
porated by four kind of bosons leading to the following 
Holstein-Primakov transformation, 



S\ = —s\ + a^a = \^y/2si — a^a 
S^ = s 2 - tfb S+ = (y/2s 2 - tfb) b 



Sq = si — e'e Sq = I y2sx — C c ] c 
S Z D = -s 2 + d)d S+ = dt ( ^2s 2 - dtd 



Introducing this representation in the Hamiltonian (^), 
and Fourier transforming to k-space, the SW Hamilto- 
nian can be written as, 



H sw — const 



/-" k V A d 



d A 



where a\. = (a k , &l k , cl k , <4j a k, &-k, c_k, dk), d = 
\diag(A 1 ,A 2 ,A u A 2 ), 
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FIG. 2. Spin wave spectrum for the 2D array. Dashed, dot- 
ted, and solid lines are for J± = 0, 0.2, and 0.5, respectively. 
For the ladder case the spectrum is qualitatively the same. 



and 



const = -4NJsls2 - N2J±(sj + s\) 
Ai = (2Js 2 + 2J ±Sl ), A 2 = (2Jsi + 2J ± s 2 ) 
2slJ ± cosky/2 ,D^ k = 2s2J±cosky/2, 
C k = 2Jy/sis 2 cos kx/2. 



We have para-unitary diagonalized pQ ] H sw which al- 
lowed us to calculate the spectrum excitation (with an 
acoustic and an optic band each one twofold degener- 
ate), the ground state energy and the magnetizations. 
For the particular case J± = wc recover the SW predic- 
tions of a ferrimagnetic chain, that is, a quadratic 
gapless dispersion near k = (Figj^) (which mani- 
fests the ferromagnetic character of the system at low- 
temperature), an antiferromagnetic gapped mode, and 
the magnetization values m\ — 0.6951 and m 2 — 0.1951 
for the spin-1 and spin-1/2, respectively. As soon as J_l 
is switched on the gapless mode becomes linear reflecting 
the antiferromagnetic character |2lJ of the coupled sys- 
tem, whereas the optical modes moves upward (Fig]2j). 
For the isotropic case J± — 1 we obtained a ground state 
energy eo = —1.3631, and mi = 0.7326, m 2 — 0.3694 for 
the magnetizations. Notice that m 2 is bigger than the 
SW value of the uniform spin-1/2 case m = 0.303. Using 
mean-field arguments, such a robust spin-1/2 magneti- 
zation is due to the more stronger field caused by the 
spin-1 neighbors (a similar reasoning explains why mi 
is smaller than the SW value of the uniform spin-1 case 
to = 0.803). It is also interesting to note that the rela- 
tive reduction of the magnetizations, m\/S and m 2 /s, are 
both nearly 0.73 in the isotropic case, while m\/S ~ 0.70 
and m 2 /S ~ 0.40 for the ferrimagnetic chain. In Fig.|| 
we show mi and m 2 versus Jj_. It can be seen that the 
effect of 2D coupling is to enhance the values of the mag- 
netizations, being the crossover completely smooth. On 
the other hand, our SW calculation can be adapted to 
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FIG. 3. Magnetizations predicted by spin wave theory. 
Solid (dashed) line represents the magnetization mi (7712) for 
spin-1 (spin-1/2) sites. 
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0.8 (empty square), 0.3 (filled circle), and 0.1 (filled square). 
Solid lines are polynomial fits as in eq.(^). 



the ladder case (a) by just changing: 



-4, 



const = 
= (2 Js 2 



-ANJsls2- NJ ± (s(^ 
-J ±Sl ), A 2 = (2J Sl - 

= slJ x ,L>4 = s2J -i- 
= 2 J^/s\s 2 cos kx/2. 



J±s 2 ) 



and k running along the ID Brillouin zone. Now, the 
linear behaviour of the gapless mode leads to important 
quantum fluctuations that prevent the ordering in ID. 
This feature was previously pointed out by Fukui et. al 
jl2| and suggests a possible opening of a spin gap which 
is impossible to predict within a SW calculation. In spite 
of such break down we have corroborated that the SW 
ground state energy compare quite well with the DMRG 
prediction. 



III. DENSITY MATRIX RENORMALIZATION 
GROUP 

The difficulty of SW theory to describe ASL leads us to 
implement the DMRG method, which has shown a great 
versatility in different quasi-onedimcnsional systems. We 
have performed calculations using both, the finite and in- 
finite algorithm, with open boundary conditions (OBC). 
The procedure follows the usual DMRG steps [^2| , think- 
ing the ladder like a chain of rungs. We found out that 
to avoid great edge effects, because of the OBC, the def- 
inition of two different single sites is needed, each one 
composed of a spin-1/2 rung and spin-1 rung, respec- 
tively (see Fig.[l](a)). The 4 and 9 states needed to span 
each rung (otherwise we would need just 6 states) ap- 
pear at the beginning as a complication to the number of 
dominant density matrix states kept, m, but they shown 
to be just of minor annoyance. Most of the calculations 
was carried out using m between 200 and 400, getting a 
truncation error O(10 -8 ) at worst and O(10 -12 ) in the 



best case, reassuring us the reliability of the calculation. 
We have calculated the triplet spin gap which is defined 
as A(L) = E{L, S z = 1) - E{L, S z = 0), with E(L, S z ) 
the ground state energy for a chain with L rungs (2 x L 
is the number of sites), and z component of total spin S z 
p3| . For each J± , we have extrapolated the gap with L 
ranging from 10 to 100 rungs using a polynomial fit of 
the form |24 |: 



A(L) = A + ai /L + a 2 /L 2 + a 3 /L 3 



(2) 



In Fig.[| it is plotted the gap A vs. Jj_. In the partic- 
ular case, Jj_ = 1, we obtained a value for the gap of 
A = 0.334. For small values of J± (below 0.3) the gap of 
the ASL behaves quadratically as A ~ 0.6JJ and beyond 
0.3 it turns out quite linear corresponding to the strong 
coupling regime. Consistently with these features, when 
extrapolating the gap for different values of Jx, we de- 
tected a crossover between two different scaling regimes 
around J± ~ 0.3. This is shown in the inset of Fig.|] 
where we have plotted the scaling for J± = 0.1,0.3,0.8 
and 1. It is interesting to compare our result for the 
weakly coupled regime with that of the USL, which is 
known to behave linearly as A ~ 0AIJ± |2j|. We think 
that this discrepancy reflects the quite distinct underly- 
ing physics behind the decoupled regime, that is, a crit- 
ical state for uniform spin chains and a broken symme- 
try state for the ferrimagnetic chains. In order to get 
an insight of the ground state configuration of the ASL 
we computed |26| the local bond strengths and compared 
them with the USL. The comparison was made after nor- 
malizing each value to its free bond case. Hence, for the 
very strong coupling regime the values of the bonds along 
the leg (rung) must go to zero (one). Actually, we cal- 
culated the z component of these quantities and because 
of rotational invariance we multiply it by a factor 3. In 
Fig.||(a), we present the normalized values of (Si.s,+i) 
along the leg. A similar monotonic behavior with J± 
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FIG. 5. Normalized bond strength as a function of Jj_ . (a) 
empty (filled) diamond along the leg for ASL (USL), and (b) 
spin-1 rung (empty circle), spin-1/2 rung (empty square) for 
ASL and spin-1/2 rung for USL (filled circle). 



is observed for ASL and USL. On the other hand, in 
Fig.||(b) we show the normalized values of (Sj.Sj+i) and 
(si.Sj+i), corresponding to the vertical spin-1 rungs and 
spin-1/2 rungs respectively. It can be noticed that the 
spin-1 rungs are always larger than the spin-1/2 ones, 
however, what is really remarkable is the crossing be- 
tween the spin-1/2 rungs values for ASL and USL which 
occurs just at J± ~ 0.3. Even if we have not a micro- 
scopic explanation for that we believe that our results for 
the gap and the local bond strengths are connected and 
indicate a change of regime for the above value of J±. 
For the sake of completeness, we studied the correlation 
functions along the legs. In Figj^(a), we compare the spin 
correlations of the USL with the three possible correla- 
tions of the mixed spin case, corresponding to J± =0.1 
p7| . In the inset it is shown in a semilog plot the expo- 
nential behaviour of the correlations with an estimated 
correlation length of £ USL ~ 25, and £ ASL — 30 lattice 
spacing. It can be noticed, also, that the three possible 
correlation lengths are the same. We have checked that 
the correlation functions along the leg decay exponen- 
tially for the whole range of interchain coupling and the 
correlation length ^ ASL — > oo in the limit J± — ► like 
in USL JIU . When J± is strictly zero there is a transi- 
tion to two ferrimagnetically ordered chains in the ASL. 
In Fig.||(b) we show the correlations for Jj_ = 1. Again, 
the semilog plot (inset) shows the exponential decaying 
of the correlations with an estimated correlation length 
of £ ASL ~ 5. This value is bigger than the isotropic one 
£ USL ~ 3 found for the uniform case I2HJ2F 
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FIG. 6. Correlation functions obtained with infinite 
DMRG algorithm, for (a) J± = 0.1 and (b) J± = 1. In- 
set: semilog plot of the absolute values. In (a) the results are 
compared with the USL. 



IV. SUMMARY AND CONCLUDING REMARKS 

In summary, we have studied the effect of AF coupling 
in mixed spin chains for two particular topologies, the 
ladder and the two dimensional system. In the latter, 
SW theory leads to a smooth dimensional crossover from 
ID to 2D with an enhancement of the magnetizations. 
Furthermore, as a consequence of the crossover, we de- 
tected an important difference in the relative magneti- 
zations, being mi/S and mij s around 0.73 for Ji = 1 
whereas mi/S ~ 0.70 and m-ij s ~ 0.40 for J± = 0. 
For the ladder system, where SW breaks down, we used 
DMRG to study the gap and the correlation functions 
with J± . For Ji 7^ the ladder system is always gapped 
with a quadratic behaviour A ~ 0.6 J\ for J± lower than 
0.3, and then it becomes quite linear for greater values 
of the interchain coupling. The study of the local bond 
strengths also indicate a change of regime at J± ~ 0.3. 
In the particular case of J± = 1 we found that the gap 
is A = 0.334 and that the correlation functions decay 
exponentially with a correlation length about £ ASL ~ 5. 
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Except for the differences found in the weakly coupled 
regime, our results suggest a strong similarity between 
the ASL and the USL. This close resemblance encour- 
age us to make some general statement regarding an ar- 
bitrary number of mixed spin chains, between the two 
limiting cases we have studied. First, we can rigorously 
say that any odd number of chains has always a ferri- 
magnetic ground state, due to Lieb-Mattis theorem, so, 
it will be gapless -but ordered-. Then, if we complement 
our results with the conjecture that an even number of 
chains will be gapped and -similarly to uniform ladders 
p8| , p4|j28| - this gap will decrease to zero in the 2D limit, it 
is recovered an analogous spin gap behaviour to the uni- 
form spin-1/2 case. In both cases, however, the nature of 
the gapless states is completely different. We hope our 
findings could be tested experimentally in a near future. 
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